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Quasimolecules, which consist of two differently excited atoms in 
a resonantly excited gas, are considered. The energy of dissocia- 
tion and typical sizes of such molecules are calculated in the first 
order of quantum-mechanical perturbation theory with the help 
of the dipole-dipole interaction operator. It is shown that there 
exist metastable quasimolecules, whose dipole radiative transition 
to the ground state (two non-excited atoms) is forbidden. The 
lifetime of such molecules is estimated and it is shown that quasi- 
molecules may considerably affect the transport processes in a 
resonantly excited gas. 



1. Introduction 

It is well known (see e.g. [I]) that the interaction en- 
ergy between atoms which do not form stable molecules 
can be calculated by means of quantum-mechanical per- 
turbation theory using a dipole-dipole interaction oper- 
ator. Particularly, the interaction between atoms in the 
ground state is described with the help of the van-der- 
Waals potential 



Uvw{R) 



const 



R:$> a, 



where i? is a distance between the atoms, and a is a typ- 
ical size of an atom. The potential UywiR) is obtained 
in the second order of perturbation theory and always 
corresponds to attraction. The above-mentioned forces 
are the so-called van-der-Waals forces, which form a rel- 
atively shallow potential well (approximately 10"'^ eV 
for He2 2 and 2 x 10~^ eV for Kr2 3J that corresponds 
to approximate temperatures 10 K and 200 K respec- 
tively), so that the formation of stable molecules in such 
gases at room temperature does not happen. 

The situation changes cardinally if one of the atoms 
is in an excited state, from which the dipole transition 
to the ground state is possible. In this case, the dipole- 
dipole interaction operator gives a non-zero contribution 
even in the first order of perturbation theory due to the 
fact that non-perturbed wave functions must account for 
an exchange of excitations between the atoms. The tran- 
sition of the excitation from one atom to the other in 



this case replaces the exchange interaction that leads to 
the chemical valence. Such an interaction between the 
atoms is conventionally called a resonance interaction. 
It is described with the help of the potential 



Ur{R) - ± 



const 



and can have a character of not only attraction, but also 
repulsion. 

The dissociation energies of bound states, which can 
be formed at the expense of the resonance dipole-dipole 
interaction, turn out to be of the order of 1 eV accord- 
ing to the below-presented estimations. This means that 
quasimolecules, which consist of atoms in the ground and 
an excited states, can be formed at room temperature. 
It should be emphasized that a relatively large energy of 
the resonance interaction is related, as mentioned above, 
to the exchange of excitations realized between identical 
atoms through the dipole-dipole interaction operator. 

It is worth noting that in the case of differently ex- 
cited identical atoms, when the dipole transitions from 
excited states to the ground state are forbidden, the 
formation of quasimolecules (so-called excimers) is also 
possible [1]. In contrast to resonantly excited quasi- 
molecules, when the attractive part of the potential en- 
ergy of the interaction as a function of R can be cal- 
culated analytically, the potential energy curves of the 
molecular states of excimer molecules, as known, can be 
calculated only by numerical methods. 

The presence of even relatively low concentrations of 
resonantly excited molecules may essentially affect the 
collective properties of excited gases and plasma. The 
thermodynamic functions of resonantly excited gases 
were first considered in [S] . Quasimolecules of two differ- 
ently excited one-electron atoms in such gases were dis- 
cussed in [5] , and a role of quasimolecules in resonantly 
excited gases was considered in [7]. 

The purpose of the present paper is an analytical 
calculation of the long-range part of the potential en- 
ergy curves of molecular states of the resonantly excited 
quasimolecules formed from identical atoms with the 
closed electron shells and atoms with one valent electron 
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that do not form stable molecules in the ground state at 
room temperature. We also show that a radiative transi- 
tion of resonantly excited quasimolecules to the ground 
state (two non-excited atoms) is forbidden in the dipolc 
approximation. Therefore, the considered quasimolecule 
is metastable and its lifetime appears to be of the order 
of 10~^ sec. 



2. Quasimolecule of Helium He2 

Consider an interaction of two helium atoms, one of 
which is in the ground state, while the other is in the 
excited P— state. Let us denote the nuclei of these atoms 
as A and B and assume that the distance between them 
equals R. We will use two coordinate systems, which 
describe the electrons of the atoms, and match their ori- 
gins with the corresponding nuclei (A and B). Next, we 
direct the z axes of these coordinate systems along the 
line that connects these nuclei, and denote the electrons 
belonging to the first atom (nucleus A) by indices 1 and 
2. The electrons belonging to the second atom (nucleus 
B) are denoted by indices 1' and 2'. 

The electron wave function of atom A (in the ground 
state) calculated with the help of the variational method 
looks like [8j 



where 



ri+r2 

*(l,2) = ^e-"^. 



(1) 



and represents the product of two hydrogen-like wave 
function for each electron in the field of the nucleus with 
an effective charge Z^s — a — 27/16 : 



Mr) = |0) = 



v3/2 



33/2 



'IT a 



3/2 



(5) 



and 



^Px (r) = Ix) = Rni (r) ^J ^ sin 6* cos ip, 
'/'y(r) = I: 



Rni{r) J ^ sine sin(/?. 



(6) 



Rni{r) cose 



The degenerated wave functions ^ are chosen such that 
they transform as vector components under rotations of 
the coordinate system. 

Henceforth, we restrict ourselves by the description 
of the interaction of atoms in the dipole-dipole approxi- 
mation. As is well known [5 , the resonance interaction 
(between atoms in the ground and an excited states) is 
inversely proportional to i?"^ and shown to be the most 
intensive for low-excited states. Therefore, we will con- 
sider the lowest excited P— state with the principal quan- 
tum number n = 2, for which 



P2i(r) = 



2 7' 



5/2 



1 



— p, 'a 'y — — 

V3 a3/2 a ^ ' ^ 2- 



(7) 



The wave function in the zero-order approximation is a 
linear combination 



*(1,2,1',2') = ^QV^,(121'2') 

i=l 

of the following 12 wave functions: 
lOOOx), 100x0), 1 0x00), 1x000), 



(2) |0002/),|002/0),|OyOO), 12/000), 



(8) 



(9) 



(10) 



Here and below, a designates Bohr's radius: 



0.529A. 



(3) 



We take the electron wave function of the excited 
P-state of a helium atom in the form of an antisym- 
metrized product of the hydrogen-like wave functions of 
an S'-electron in the field of a nucleus with charge Z = 2 
and a P-electron in the field of the atomic residue with 
charge Z = 1 [S]: 



VI/(1',2') = i= [M^')^^{'2')±M^>^{^')] , ^=X.y,Z, 



lOOOz), lOOzO), lOzOO), |zOOO). (11) 

Wave functions (P|)- pT|) arc the products of the one- 
particle wave functions ([2]), ([5]), and ([6]) with the electron 
coordinates taken in the order 121'2', e.g. V6(121'2') = 

lOOyO) = Mri)Mr2)Vy{Yv)Mr2'). 

The interaction operator in the dipole-dipole approx- 
imation has the form 



V 



(di -hd2)(di/ +d20 - 3(di +d2,n)(di/ +d2',n) 



i?3 



(12) 



where n = R/P, and the dipole moments of electrons in 
their proper coordinate systems are equal to 



(4) di, 



-eri,2, dy 



(13) 
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In the first order of the stationary perturbation the- 
ory the secular equation 



detl^y - AES,j\ = 



(14) 



(determinant of the 12-th order) constructed of matrix 
elements of the interaction operator on functions ([9])- 
(fTT|) brakes into 3 blocks for each of the set of functions 
(l9])-(IIl|). Namely, for sets ©-(Uni) we have 



-AE 





A 


A 





-AE 


A 


A 


A 


A 


-AE 





A 


A 





-AE 



= 0, 



and, for set (|lip . we have a similar equation with value 
—2A in the place of A. 

Here, constant A corresponds to 24 non-zero matrix ele- 
ments of the dipole-dipole interaction operator on wave 
functions (P|)- (fTT|) . This value can be calculated as 

A^ (60a;0|y|2:000) ^ ^ {0\0)^ (a;|d^|0)2 
-ft'' 



e 



dV 



(a/3) 



3/2 „-(a+/3); 



-I 2 



dV 



v3/2, 



216 Q,6/53^5 



tfx {r, 9,ip) r sin 9 cos (p 



i?3 (a + /3)6(a + 7) 
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The substitution of constants a 
7=1/2 gives 

^2„2 

A ~ 0.06- 



27/16, (} 



i?3 



(16) 
2, and 

(17) 



From the secular equation, which reduces to 3 equa- 
tions of the 4-th order, namely to AE^ - AA^AE"^ = 
for x, y-functions ©, ^ and to AE^ - IQA^AE'^ = 
for 2-functions pT|) . we obtain the energy of states, into 
which the 12-fold degenerated energy level of the system 
of two non-interactive atoms is split: 



AE = ±2A,0,0 
for X, ^/-functions and 
AE = ±4A,0,0 
for z-functions. 



(18) 



(19) 



The 6-fold degenerated level with zero energy and the 
levels with positive energies (-I-2A, -I-4A) are "dissocia- 
ble" levels and they do not form any bound states. We 
are interested in the levels with negative energies: the 
doubly degenerated (in x and y) level with the energy 
—2A and the wave function (a linear combination whose 
coefficients are eigenvectors of the correspondent secular 
equation) 



(x) 



1 



2A = 2 ( 1000^) 



|00a;0) - lOxOO) - 1x000)) , (20) 



(15) *-'2A - ^ ( lOOOy) + lOOyO) - lOyOO) - |yOOO)j , (21) 



1 / 



and the non-degenerated level with the energy —AA and 
the wave function 



(2) 

-iA 



1 / 



( lOOOz) + lOOzO) + iGzOO) + |zOOO)j . (22) 



The eigenvalues (energies) are written as indices in wave 
functions (|20 l) - (p2)) . 

According to the classification of electron levels of 
diatomic molecules (see e.g. [1]), the considered excited 
P-states of a quasimolecule with energies —2A and -A A 
belong to ^Uu and ^Eg states, respectively. 

Let us now consider determination of the equilibrium 
distance R between the nuclei. It is impossible to find it 
by accounting for only the long-range dipole-dipole at- 
tractive forces. Therefore, we have to consider repulsion 
forces that act at small distances between the atoms. 
We calculate the equilibrium distance using the follow- 
ing simple model: we assume that the electron shells 
of two atoms are not deformed by the interaction, and 
that the equilibrium distance is attained only when the 
maxima of their electronic probability densities are in 
contact with each other (notice that even this configura- 
tion corresponds to the essential overlapping of the wave 
functions). 

The maximum of the electron probability density of 
a helium atom in the ground state can be found from 
the equation 

" '^^ 0, 



d_ 

dr 



d 
dr 



a 



2a- 



which gives 

^max = ci/ct = 16/27a ~ 0.6 a. 



(23) 



(24) 



In a similar way, we find the maximum of the elec- 
tron probability density of a helium atom in the ex- 
cited P— state. In this case, the most probable value 
of r equals 



4a, 



(25) 
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but the wave function, in addition, depends on the angle 
coordinates 0, if. Accounting for that, we find the max- 
imum of the distribution of the the probabihty density 
of a;, j/-wave functions Yp\\^ = \ fx\'^ + \ fy\'^ to be at 



(26) 



''max — r* sin 9 = Aa sin 



From contact condition of the surfaces (p4)l and ([26]), we 
find the unknown equilibrium distance: 



R^''-'^^ ~ 1.203 a, 



28.4" 



Similarly, for z-functions wc have: 
Tmax = r* cos^ 9 = Aa cos^ 9, 
wherefrom we get 
i?^^^~4.6a, (61 = 0°). 



(27) 



(28) 



(29) 



At last, knowing the equilibrium distances between 
the nuclei, we are able to calculate the corresponding 
binding energies (dissociation energies with the opposite 
sign): 



AE^^'^y^ = -1.875 eV, 



-0.067 eV. 



(30) 



(31) 



We must emphasize that the total energy of a quasi- 
molecule relative to the energy of non-interacting helium 
atoms in the ground state is larger than the calculated 
dissociation energy by the value (experimental data j9j): 

E{ls2p)^E{ls^) ~ -57.787 eV-K 79.005 eV ~ 21.22 eV. 

(32) 

This energy must be used in considerations of radiation 
energy of quasimolecules (see the next section). 

3. Lifetime of Hej quasimolecules 

Let us now discuss the lifetime of quasimolecules. For 
such a study, we have to use the radiation theory of 
molecular systems. The total probability of radiation 
is written in the form of a perturbation theory series, 
each member of which represents a multipole radia- 
tion of a certain type (electric or magnetic) and order 
(dipole, quadrupole etc). Moreover, it is well known 
that the increase of the order of multipolarity by 1 
decreases the radiation probability by ~ (fcoo)^ times, 
where ag is a characteristic size of the radiative system 



and fc = 27r/A = w/c is a wave number of the radiated 
light [10,. 

In our case of a quasimolecule in the state ^n„, the 
value of kao is of the order of 10~^ {huj ~ 20 eV, 
flo ~ lA), which gives the decrease of the radiation prob- 
ability in 10** times when we proceed to the radiation 
with multipolarity by 1 order higher. For example, a 
characteristic time of the dipole radiation of an excited 
atom is of the order of 10~^ sec, and in the case where 
a dipole transition is forbidden, the lifetime increases by 
4 orders, i.e. up to 10"^ sec. Let us show that such a 
situation indeed takes place for helium quasimolecules 
Hcj in the state ^n„. 

A dipole matrix element for the transition of a quasi- 
molecule to the ground state (two free atoms) is 

dfi = (*(-^^^^)| di+d2+dp+d2HV^o(l)^o(2)^o(l')^o(2')), 

(33) 

where vpl^.f ^^i) ^^le correct wave functions of the zero- 
order approximation (pUjl - ip^ . 

It is easy to see that after the interchange of electron 
pairs 1,2 and l',2' which describe physically identical 
states, wave functions (1^01) and (PT|) change their signs 
as well as the integral As a consequence, we get 

that this integral is equal to zero 

(Mr(-,?^)| di+d2 + di,+d2.|V^o(l)^^o(2)^^o(l')^o(2')) -0, 

(34) 

and, therefore, the dipole transition from the state 
to the ground state is forbidden. 

However, for the state (z-functions ((22)l ). the 
symmetry does not forbid the dipole transition. In this 
case, only the matrix element of the z-component of the 
dipole moment is non-zero and turns out to be 

(*(-)| (di +d2+ dy + d2').|V^o(l)^o(2)V^o(l')^o(2')) 
29 „3 ^3/2 ^5/2 



■ e a ~ —0.49 ea. 



(35) 



The probability of a dipole radiation of the excited sys- 
tem per unit time equals 



(36) 



wherefrom one can estimate the lifetime of the state: 



1 3rc' 



4^3 



w 4(fiw)M0-49ea)2 



0.18 ns. 



(37) 
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We see that the ^Ylg state is short-hved and cannot a 
metastable one in contrast to the ^n„, whose dipole ra- 
diation is forbidden and therefore a quasimolecule in this 
state is metastable and lives at least 10~^ sec, that is es- 
sentially greater than characteristic times of the dipole 
radiation 10^® ~ 10~^ sec. 

Let us note, that the given estimation (|37p of the 
quasimolecule lifetime agrees (by the order) with the ex- 
perimental data for short-lived excited states of helium 



of the zero-order approximation of two non-interactive 
atoms with account for the excitation exchange are con- 
structed similarly to what was done in Section!^ Namelv. 
they are 

\0x), \xO), |0y), |yO), |0z), |zO), (39) 
where 

Mr) = \0)^Rno{r) (40) 
V47r 



4. Quasimolecules of Atoms with One Valent 
Electron 

Let us consider a system, which consists of two identical 
atoms, and where each atom has only one outer electron 
while all the other electrons form entirely occupied elec- 
tronic shell (chemical elements of the first group). We 
assume that one of the atoms is in the ground S'-state 
(i.e. the valent electron has ns^-configuration), but the 
other one is in the excited P-statc (np^ -configuration). 
The total electron wave function of each atom can be 
represented in the form of the product of the electron 
wave function ^'(ri, . . . , rjv) of entirely occupied elec- 
tronic shell and the wave function ?/;s^p(r) of the outer 
electron: 



* = *(ri,...,rAr)V's,p(r), 



(38) 



where = 2, 8, 10, 18, . . . is the number of electrons in 
the entirely occupied electronic shells. 

The entirely occupied electronic shells have zero to- 
tal spin and zero total angular momentum. Therefore, 
they do not take part in the chemical combinations and 
the chemical valence is conditioned only by the outer 
electrons [T]. Hence, the wave function ^'(ri, . . . , rjv) 
appears invariantly in all the equations where chemical 
combinations are considered and practically drops out 
from calculations. Taking this fact into account, we will 
further use only the wave function of the valent electron 
V's,p(r), but will keep in mind that the entirely occu- 
pied electronic shells must be considered for obtaining 
equilibrium distances between the nuclei of the atoms. 

Thus, the problem reduces to the interaction of 
two hydrogen-like atoms, where one electron moves in 
the field of an atomic residue with effective charge 
Z{r), which can be calculated if the wave function 
5'(ri, . . . , Tat) is known. But such a procedure can be 
done only by numerical methods. 

The angular part of the wave function of the outer 
electron is the spherical function Yira{0,Lp), but the ra- 
dial one Rni{r) is considered unknown. Wave functions 



V'p(r) 



Rn'i{r) y 3^ sine cos(p, 
Rn'i{r)\[^ sin 6* sin.^, 

J3_ 

4;t 



n' > 2. 



\z)=Rn'i{r)J^ COS0. 



(41) 



We emphasize that we neglect the exchange interaction 
of the outer electron with the inner ones. 

The dipole-dipole interaction operator is written in 
the form 



V = 



(didi0-3(di,n)(di>,n) 
i?3 



(42) 



where indices 1,1' correspond to electrons that belong 
to the different atoms with the distance between them 
being equal to R. 

After solving the secular equation constructed with 
the help of interaction operator (|42|) and the set of de- 
generate wave functions (f39|) . we get 



where A is 

^2 



^ " 3^ (/ Rnoir)Rn'iir) dr 



(43) 



(44) 



In the case of a hydrogen atom, the radial functions 
Rni{r) are known. Particularly, for n — I and n' ~ 2, 
we have 



A ~ 0.555 



i?3 ■ 



(45) 



With the help of the method described in Section [23 we 
estimate the equilibrium distance between the nuclei 

jl{^,y) ^ 1J6 a (61 « 32.8°), = 5 a {0 = 0°), 

and get the following binding energies of quasimolecules 

^ _2.77eV, ~ -0.24eV (46) 
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in the states and '^Sg, respectively. 

For other atoms, the radial functions are unknown, 
and we use the experimental atomic radii rat for es- 
timations of the dissociation energies of quasimolecules 

2 2 

and assume that A « The results of calculations 

(by the method of Section [2]) for n = n' are listed in the 
Table: 



Atom 


n 


rat, A 


a 




a 


E'-^>, oV 


Li 


2 


1.520 


3.70 


-0.54 


5.75 


-0.29 


Na 


3 


1.858 


4.53 


-0.29 


7.02 


-0.16 


K 


4 


2.272 


5.54 


-0.16 


8.59 


-0.09 


Rb 


5 


2.475 


6.03 


-0.12 


9.36 


-0.07 


Cs 


6 


2.655 


6.47 


-0.10 


10.04 


-0.05 



The typical sizes of quasimolecules can be roughly es- 
timated by the formula Rm — '2R^^\ Hence, the typical 
sizes of quasimolecules from Li to Cs lie in the interval 
from 6 to 10 A. This means that the corresponding ef- 
fective scattering cross-sections are approximately in one 
order of magnitude larger than the gas-kinetic ones. 

As well as for helium quasimolecules, the dipole tran- 
sition appears to be forbidden in the case of the atoms 
with one valent electron due to the fact that the symme- 
try of the wave function remains unchanged with respect 
to the exchange of the nuclei. Thus, the quasimolecules 
are metastable only in the a;, y— states. 

5. Conclusions 

The analytic calculations of the molecular states of di- 
atomic quasimolecules, which are formed of the atoms in 
the ground and an excited states, are carried out. The 
calculations are performed in the first order of pertur- 
bation theory with the use of the dipole-dipole interac- 
tion operator for atoms with entirely occupied electronic 
shells and atoms with one valent electron. The results 
describe molecular states at the distances R greater than 
the typical size of atoms. At distances (phe- 
nomenological parameters of the theory) of the order of 
atom sizes, where the repulsion of atoms predominates, 
the potential energy is modelled by a solid wall. The 
dissociation energy of such molecules, according to our 
estimation, appears to be of the order of 1 eV. In con- 
trast to excimer molecules, whose typical lifetime varies 
from 0.1 to 10 ns [4], quasimolecules are metastable. The 
estimation shows that their radiation decay is forbidden 
for dipole transitions and is of the order of 10~^ sec. 

The above-described quasimolecules exist in active 
media that are used for the generation of laser radiation 
[1], [2] J [II]- It is worth noting that, even for the sim- 
plest excimer molecule He2 there exist over 60 molecu- 
lar potential curves [1] . This fact essentially complicates 
theoretical analysis of systems containing a noticeable 



quantity of excited atoms from the point of view of iden- 
tification of certain spectral levels. 

The presence of excimers and quasimolecules essen- 
tially affects collective properties of the excited gases. 
For example, one can affect the diffusion and the thermal 
conductivity coefficients in excited gases with regard for 
the fact that quasimolecules have much bigger effective 
cross sections than the gas-kinetic ones. In particular, 
the presence of quasimolecules leads to a deceleration of 
the processes of diffusion and thermal conductivity. 

The presence of diatomic quasimolecules may also as- 
sist a more active process of formation of clusters which 
consist of many atoms [13 . The formation of quasi- 
molecules with a finite lifetime r and their decay (disso- 
ciation) can be considered as a photo-stimulated chem- 
ical reaction. As known [14j . the anomalous sound ab- 
sorption may be observed in this case at the expense 
of the mechanism of the second viscosity at frequencies 
LUT ~ 1. Furthermore, a theory of non-equilibrium phase 
transitions in many-particle systems with finite lifetime 
is being intensively developed T,5j. The most essential 
conclusion of such a theory is that there exist additional 
restrictions on the decay of systems into two coexistent 
phases (in our case — phases with different concentra- 
tions of excited atoms), which requires additional im- 
provements [16] . 
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